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Abstract. 

We consider the post-Newtonian Umit of a general class of bimetric theories of 
gravity, in which both metrics are dynamical. The established parameterised post- 
Newtonian approach is followed as closely as possible, although new potentials are 
found that do not exist within the standard frame- work. It is found that these theories 
can evade solar system tests of post-Newtonian gravity remarkably well. We show 
that perturbations about Minkowski space in these theories contain both massless 
and massive degrees of freedom, and that in general there are two different types of 
massive mode, each with a different mass parameter. If both of these masses are 
sufficiently large then the predictions of the most general class of theories we consider 
are indistinguishable from those of general relativity, up to post-Newtonian order in 
a weak field, low velocity expansion. In the limit that the massive modes become 
massless, we find that these general theories do not exhibit a van Dam-Veltman- 
Zakharov-like discontinuity in their 7 parameter, although there are discontinuities 
in other post-Newtonian parameters as the massless limit is approached. This smooth 
behaviour in 7 is due to the discontinuities from each of the two different massive 
modes cancelling each other out. Such cancellations cannot occur in special cases with 
only one massive mode, such as the Isham-Salam-Strathdee theory. 



t JS Guggenheim Memorial Foundation Fellow 



The Parameterised Post- Newtonian Limit of Bimetric Theories of Gravity 



2 



1. Introduction 

There exists a long liistory in gravitational physics of considering theories that are 
generahsations of, or alternatives to, general relativity. These theories can take on a 
variety of different forms, and often involve the introduction of extra scalar [1] , vector |2] 
or tensor fields [3] in the gravitational action. Here we consider the latter of these: 
Bimetric theories in which there are two dynamical rank-2 symmetric tensor fields, rather 
than the usual one. Our goal is to calculate the Parameterised Post-Newtonian (PPN) 
limit of a general class of these theories, so that their consequences for observational 
and experimental gravitational physics can be determined in a straightforward way. 

The PPN formalism [U |5] is now the standard frame-work within which 
investigations of the phenomenology of relativistic gravitational physics are performed 
[6]. It has at its heart a simple and compelling rationale: that by observationally 
constraining terms in a general 'test metric', it is then possible to impose corresponding 
constraints on a variety of different metric based gravitational theories without going 
through the rigmarole of directly calculating their physical consequences every time. 
This allows the space-time geometry of the physical environments in question to be 
constrained in a theory independent way, as well as supplying a quick and direct route 
to determining the observational validity (or otherwise) of particular theories. 

The problem of determining the PPN limit of bimetric theories of gravity has been 
approached in the past. The PPN limit of the bimetric theory of Rosen [7] was found 
by Lee, Caves, Ni and Will [8]. The PPN limits of Rastall's bimetric theory piTO]. 
and Lightman and Lee's bimetric theory [11] are also known. These theories, however, 
contain only one dynamical rank-2 tensor each, the other being a priori specified as 
being Riemann fiat. Such non-dynamical fields can cause serious problems for these 
theories, particularly with respect to the emission of gravitational waves from binary 
systems. For example, Rosen's theory predicts an unobserved increase in the rotational 
period of the Hulse- Taylor binary pulsar system PSR 1913-1-16 [12], and it can be argued 
that similar behaviour should be expected in all theories with 'priori geometry' [6] . The 
bigravity theories considered in the present article do not have any such non-dynamical 
fields, and are therefore not expected to fall foul of binary pulsar observations in the 
way discussed above. To the best of our knowledge, the case of determining the PPN 
limit of a general class of bimetric theories with two dynamical rank-2 tensor fields has 
yet to be investigated in detail. It is this subject that we intend to investigate here. 

As motivation for this study we use the recent astronomical observations that have 
revealed that most of the matter in the Universe is in the form of unknown "dark 
components" (i.e. particles or fields that are not present in the standard model, and 
that do not interact electromagnetically) . Despite ongoing efforts these fields have yet 
to be observed directly, and it seems natural to consider extensions of general relativity 
that may explain them. In this context, bimetric gravity has recently been put forward 
as an alternative to both dark matter and dark energy. For details of the way in which 
the extra degrees of freedom in these theories can be made to reproduce the observations 
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usually attributed to the dark components of the Universe we refer the reader to [T3HIS]- 
Here we will not concern ourselves with the cosmological implications of these theories 
directly, but will instead study their weak field limit. If theories of this type are to be 
considered observationally viable it is absolutely necessary that we properly understand 
their weak field limit, and how they are constrained by observations of post-Newtonian 
gravitational phenomena. 

The plan of this paper is as follows. In Section |2] we introduce and discuss the 
bimetric theories of gravity that we will be investigating, giving their gravitational 
action and field equations. In Section [3] we briefly introduce the PPN formalism, the 
perturbative expansion it relies upon, and the way in which it is used to constrain 
gravitational theories. Section H] contains a calculation of metric perturbations up to 
Newtonian accuracy. We then proceed in Sections [5] and [6] to derive the full post- 
Newtonian limit of the metric perturbations about both of the metrics involved in 
these theories. In Section [7] we transform into the standard post-Newtonian gauge, and 
calculate the form of the perturbations in the metric that couples to matter. Section [8] 
provides a discussion of the extent to which these theories can be subject to observational 
constraint, and in Section |9] we provide a brief discussion of our results. 

2. Bimetric Theories of Gravity 

One of the most straightforward ways to extend general relativity is to take the already 
existing concept of a dynamical rank-2 symmetric tensor field, and replace it with a 
multiplet of Ng such fields. Multigravity theories can then be constructed by considering 
an extension of the Yang-Mills approach to gauge theories. That is, if we write our Ng 
fields as g'^^^, where a = 1, 2, A^^, then our gravitational action can be written 



where lint is an interaction term that can depend on any or all of the metrics. Although 
it is impossible to construct a non-trivial lint term that preserves the original A^^- 
dimensional symmetries of the free theory [20], the action above is manifestly invariant 
under diagonal diffeomorphisms acting on all metrics. The potential therefore breaks the 
symmetry group down to the diagonal subgroup of diffeomorphisms acting on all metrics 
with the same parameter. This approach to gravity has a long history, especially the 
Ng = 2 incarnation known as 'bigravity'. It began in the early 1970s with the pioneering 
work of Isham, Salam and Strathdee [3] , while today actions like ([T]) are mostly used as 
covariant actions for massive gravitons. A classification of different interaction terms, 
as well as a variety of different motivations for bigravity, are explored in [2T] . 

In order to make progress we must specialise the general action ([T]) somewhat. The 
theory we will consider for the rest of this paper is therefore the Ng = 2 bigravity theory 




(1) 
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^(i?-2Ao) + v^(ir-2A) 



P 



-q [koJ + KiJ^ + K2{q-T{q-')^,) 



(2) 



Here we have written the two metric fields as gfj,^ and qf^u- The Ricci scalars constructed 
from these metrics are written as R and K, respectively, and the two constants Aq 
and A are the bare cosmological constants on g^^, and q^^,. The three constants Kq, 
Ki and K2 parameterise the interactions between ^f^^ and q^^. The inverse {q~^)^'^ is 
defined such that {q~^Y"q^i3 = S"^, and all raising or lowering of the indices of these 
tensors is done with g^^, so that e.g. {q~^)^u = {.(l~^)°''^gafigi3u- We have also defined 
J = {q-^Y'^g^, = {q-^Y^. 

The theory above is a generalization of the often considered Isham-Salam-Strathdee 
(ISS) action 

^ ^ " ^ ^{R - 2Ao) + ay^ik - 2A) - ^y^hnt 



S 



d X 



(3) 



i9,.-f,.)i9a^-faMf-'Y''if-'r^-if-'rif-'r^)- 



16nG 

with the interaction term I^nt 
Here the scalar K is the Ricci scalar constructed from the metric f^^, and {f^^Y^ffj.p — 

The generalised theory ([2]) can be seen to reduce to the ISS action ([3]) under the 
metric rescaling q^^, = af^y, together with the identifications = 6, ki = —a and 
1^2 = 0", and the redefinition A = A/o" — 6/(Z^(j). 

Of course, we also want to couple matter to these theories. We do this by including 
a Lagrangian density of the form = Cm{ip, g^J.u), where the matter fields ip are coupled 
to a linear combination of the metrics 

9f,u = mg^y + ng^^, (4) 

where m and n are constants. After varying with respect to each of the metric fields 
this then results in the field equations 




(5) 



(ko + 2kiJ) (^q ^)^y - ]j^g^lu■J^ 



+ 2k2 ( {q-'),M-')\ - ^{q''U{q-'T'9,. 



and 



fK^, - SirGl'nJ^ ( q^aquf^T'^" - \q^.uT''^qap ) - m^.. (6) 



^ [kiJ^ + K2{q-^)ap{q~^T^] q^,u + N + 2ki J] + 2/t2(g"^; 
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where = AqP and a = XP, and where T^'^ is defined by 

Expressing the field equations in this form, without R and K, the perturbation equations 
take on their simplest form. 

Now, we should be aware that these theories contain within their spectrum a 
massive spin-2 field, and that such terms often have associated with them certain 
discontinuities and instabilities [22] - [28] . These include the van Dam-Veltman-Zakharov 
(vDVZ) discontinuity [29l[30], and the Boulware-Deser (BD) instability [31j. Both of 
these are potentially very serious problems. The vDVZ discontinuity involves the zero- 
mass limit of massive theories of gravity. Instead of recovering the solutions of general 
relativity in this limit, it can instead sometimes be found that space-time geometry 
approaches a limit which is inconsistent with observed gravitational phenomena. For 
example, the prediction for weak lensing by the Sun in the Pauli-Fierz theory [32] yields 
results which disagree with general relativity by 25%. The BD instability involves 
the possible existence of negative energy states at the non-linear level, and a lack of 
boundedness from below. We refer the reader to [331435] for recent discussions on massive 
(single metric) gravitons. 

While deserving of careful study, these problems are not necessarily fatal for the 
bimetric theories we are considering. The vDVZ discontinuity and the BD instability 
were both identified in the context of bimetric theories in which only one of the metrics is 
dynamical (the other being a priori specified as Riemann fiat). It remains to be seen to 
what extent these problems endure when both tensor fields are dynamical. Furthermore, 
in the case of the vDVZ discontinuity it is thought that in the zero-mass limit the 
graviton can become strongly-coupled, such that the usual perturbative analysis can 
break down on scales smaller than the Vainshtein radius, so that general relativity is 
recovered [36l|37]. This 'Vainshtein mechanism' can be demonstrated numerically [38] . 
Also, in the theories we are considering the vDVZ discontinuity is not necessarily a 
problem as both massless and massive modes can be shown to exist when Minkowski 
space is perturbed. The mass parameters can therefore be assumed to be large, and there 
will still always exist long-ranged massless modes to carry the Newtonian gravitational 
force (as long as matter doesn't couple to the massive modes only). The relevant question 
in this case then becomes whether or not general relativity can be recovered in the limit 
M — )■ oo, rather than M — )■ 0. This is primarily the question we will concern ourselves 
with here, although we will also pause in the final sections to consider the zero-mass 
limit of the theory. 



3. The PPN Approach 



Here we will briefiy recap the essential elements of the PPN formalism, as required for 
the internal coherence of this article. The PPN formalism is a perturbative treatment. 
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and requires a small parameter to expand in. An "order of smallness" is therefore 
defined by 

f/~t;2~^~n~ 0(2) (8) 

where U is the Newtonian potential, v is the velocity of a fluid element, p is the fluid 
pressure, p is its rest-mass density and 11 is the ratio of energy to rest-mass densities. 
We also take time derivatives to add an extra order of smallness each: 

|~0(1). (9) 

The PPN formalism then proceeds using an expansion in this order of smallness. For 
our theories, we will now perturb our two fundamental tensors as 

g^Lu = Vf^u + (10) 
Qf^y = f]f,u + hf,u, (11) 

where rj^^ is the metric of Minkowski space, and where we take f]^,^ = XqTi^^, where 
Xq is a constant. This last equation expresses the fact that, in general, one set of 
coordinates cannot be used to correspond to the same proper separation between events 
in the geometries associated with each of the two metrics. To be fully general one could 
consider different constants in front of the different components of rjn^. We will not do 
this here. 

Let us note that Minkowski space is not always a solution to the field equations ([5]) 
and ([6]), due to constant terms in the action. For q^^ = X^g^y these terms correspond 
to the 'cosmological constants' terms 

Pa = (ao + X^/to + 8^1 + 2^2) (12) 

1 f : f^o 



for g and q, respectively. These can, in general, be non-zero even in the absence of a 
and tto, resulting in a de Sitter backgrouncQ. For vanishingly small pA and px, however, 
we can approximate this as Minkowski space. This corresponds to the conditions 

ao - X^Ko - 8ki - 2k,2 (14) 
a ^ - (15) 

which we will use to eliminate ao and a in the perturbed equations that follow. 

Now consider the energy- momentum tensor, T^'^. For a perfect fluid this is given 
by T'^'^ = [p{l + n) + p] u^u'^ +pg^^ , where p is the density of rest mass, p is pressure, 11 
is internal energy per unit rest mass, and is the 4- velocity of matter. Neglecting 0(1) 
contributions from hoi and hoi, the time-like normalisation u^u'^g^^ = — 1 then gives 



An exception is the ISS theory, for which pA = ao/'&T^Gt 



•1 
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to the relevant order in perturbations the components of the energy-momentum tensor 
are then given by 

(m/ioo + ^^oo) 



rpij 



P 

[m + nXl) 
[m + nXl) 



1 + n 



(m + nX^] 



+ 0(6) 



+ 0(5) 
+ 0(6). 



(17) 
(18) 
(19) 



(m + uXq) 

The two expanded metrics f fTOl) and ffTTl) can now be substituted into the field equations 
and along with f[T7|) . ffTSl) and ffTIJl) . The field equations can then be solved for 
order by order in smallness of perturbations, and gauge transformations of the form 

^x^' + (20) 

can be used to transform it into the "standard post-Newtonian gauge", where the 
spatial part of the metric is diagonal, and where terms containing time derivatives 
are removed wherever possible. The metrics that result then allow us to determine the 
post-Newtonian parameters of these theories by comparing to the test metric 

-1 + 2GU - 2(3G^U^ - 2iG^^w + (27 + 2 + + Ci - '^i)G^i 
+2(1 + 37 - 2/3 + C2 + 0G"$2 + 2(1 + C3)G<I>3 
+2(37 + 3C4 - 20O<l>4 - (Ci - 2i)GA 



900 



-^(3 + 47 + ai - a2 + Ci - - ^(1 + ^2 - Ci + 20GW, 



(21) 
(22) 
(23) 



9oi 

g,, = il + 2^GU)6,,, 

where (3, 7, ^, (i, (2, Cs? cn, ^2 and 0:3 are the post-Newtonian parameters, U is 
the Newtonian gravitational potential that solves Poisson's equation, and $vf, ^i, $2, 
$3, $4, A, Vi and Wi are the post-Newtonian gravitational potentials given in [6]. The 
particular combination of parameters before each of the potentials in ( 12T|) . ( l22il and ( l23ll 
are given so that the parameters have particular physical significance, once gravitational 
phenomena have been computed. We have chosen to construct the test metric here in 
terms of the combined metric g^^, as this is the metric that couples to the matter fields, 
and hence is the one being constrained by observations of those fields. 

In general relativity the PPN parameters in (l2Tl)-(l23l) are given by /3 = 7 = 1, 
with all other parameters equaling zero. The interpretation that is often given to these 
parameters is that 7 is 'the spatial curvature per unit rest mass', /3 is 'the degree of 
nonlinearity in the law of gravity', the are due to 'preferred- frame effects', and ^ 
is due to 'preferred-location' effects. The parameters Q, as well as 03, are sometimes 
associated with the violation of conservation of momentum. These parameters are all 
constrained by observations, and for the present study the two constraints below are of 
particular interest [39] : 

7- 1 = (2.1 + 2.3) X 10"^ (24) 
\a2\ < 1.2x10-^. (25) 
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The first of tliese comes from observations of the Shapiro time delay of radio signals from 
the Cassini space-craft as it passed behind the Sun [39], and the second is derived from 
the observation that the Sun's spin axis is closely aligned with the angular momentum 
vector of the solar system |?T]. For further constraints the reader is referred to [6]. 

Finally, the equations of motion show that for time-like particles propagating along 
geodesies the Newtonian limit is given by c/qq (and hence goo and goo) to 0(2) and c/oi 
to 0(1), with no other knowledge of the metric components beyond the background 
level being necessarjjJJ]. The post-Newtonian limit for time-like trajectories requires a 
knowledge of 



5'oo 


to 


0(4) 


(26) 


goi 


to 


0(3) 


(27) 


9ij 


to 


0(2), 


(28) 



where Latin letters are used to denote spatial indices. To obtain the Newtonian limit of 
trajectories followed by null particles we only need to know the metric to background 
order. The post-Newtonian limit of null trajectories requires ^^oo and gij both to 0(2), 
as well as goi to 0(1). 

In order to proceed in calculating the form of the weak field metric g^^, to the 
orders of smallness specified above, it is useful to have expressions for various perturbed 
quantities involving the metric tensors, such as their determinants and inverses. These 
are given in Appendix A. 



4. Newtonian Perturbations to 0(1) and 0{2) 



4-.1. The goi and goi terms, to 0(1) 

Unlike the case of gravitational theories with a single metric, for the theories being 
considered here we cannot automatically set both goi = and qoi = to 0(1) through 
gauge transformations. We must therefore calculate these terms explicitly. 

Substituting the expressions derived above into the — i components of the field 
equations ([5]) and IQ, and eliminating ao and a using f|T4|) and f|T5l) . gives 



I' 



Oi 



2 



8Ki AK2 

8^1 4/€2 
^0 + ^ + ^ 



[Xq/ioj — hoi\ ■ 



(29) 
(30) 



The left-hand side of these equations vanishes when acted upon with V-, due to the 
Bianchi identities. We must therefore have hoi,i = X^hoi^i- The equations (l29l) and (!30ll 
can then be decoupled by defining two new variables: 



II The 0(1) terms in goi are usually omitted from the beginning in theories with a single metric, as they 
can be removed by a suitable gauge choice. In what follows we will find that the situation is somewhat 
more complicated in multigravity. 
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to get 



- ho.. - ^ (31) 
=hoi + hou (32) 



V^/iJ,? = M'^hP (33) 



V=C = (34) 



where we have defined 



M' ^ ^^^^ f°^ [Xo^«:o + 8«:i + 4^:2] • (35) 

We can now set = by a gauge transformation; the h^^^ mode, however, is gauge 
invariant, and so cannot be made to vanish in this way. We now note that h^Q^j = 0, so 
/ig™'* is the curl of some vector potential Ai, with the solution 

-M\x-x'\ M\x-x'\ 

= (V X A), = (cO.^ J- + {c,).j jr. (36) 



UM en then is a hyperbolic function, and if M G X then it is oscillatory. This 
solution does not depend on the matter content of the space-time, and exists in a vacuum 
where only gf^i, and g^jy are present. It also has a preferred point in space, x', and 
being the curl of a vector, looks like a rotation field. Finally, we impose C2 = as a 
boundary condition, in order to maintain asymptotic fiatness at infinity. This is one of 
the PPN assumptions, up to cosmological terms. 

4.2. The g^oo and goo terms, to 0(2) 

We will now proceed to calculate goo and goo to 0(2). To do this we will need both the 
0-0 field equations to 0(2), and the trace of the i-j equations to 0(2). As before, in 
order to decouple the equations let us now define massive and massless combinations 

- V - (37) 
^0 

/iW =h^, + h^,. (38) 

With the gauge choice /iq°^ = 0, and by using h^j = 0, we can then write the relevant 
massless combination of (|5]) and (|6]) as 

V^^? + 87rG(m + nX^fp 



(1 + X2) V 0*'^' °^'V 2(1 + X2 

while the massive combination is given by 
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2/,(m) 



00 



M^/ij,? + 8nG{m + nX^){m -n)p 



, (m) , (m) 
"■OiJ "'Oj,i 



Am)2 



+ 



(Xo%-8(l-X2)/s:i 



+ 
2(1 



2(2/s:i + /t2)(l+Xg 

"Oi ■ 



-A 



Here we have defined 



^ — "a "-00 ' 



(40) 



(41) 



which is given by the solution of 

{xIkq + 12ki + 6K2) - noM^XlA 



SnGXlMH^im + nXl){m - n) AM'^{8ki + (3 - X2)k2 



P + 



+ 



+ 



(X2(l - X^)ko + 8(3 - X2)«:i + 8(1 - X^)k2) , („)2 



(^0^(1 



(1+^0^) 

X2)«:o + 8(1 -X2)«:i + 8/^2) 



1 + ^0') 



(m)2 
Oi 



, (m) , (m) 



(1 + ^0^) 

To find the equations above we have made use of the expression 
^"^^'^ as well as the Bianchi identities to 0(2): 

(„) _ [X2(l - Xi)^o + 8(1 - X2)/€l + 2(3 - X2)/€2 



(42) 



h 



24-f + 2M2/io, , 
MH^Xl 



(m)2 
''Oi 



2(1 + X2; 



h 



(1 + ^0') 

[X2(l - Xi)Ko - 8XiK, + 2(1 - Xi)K2] 



Am) Am) 
"'Oj,i"'Oi,j 



2(1 + X2) 
+ l(Xo\o + 4«:i + 2/s:2)V^ [h^^ - h 



(m)2 
*0i 



(m) 
00 



(43) 



It is interesting to note that in the ISS case the factor multiplying the differential 
operator in P2|) vanishes. The equation above then becomes an algebraic relation 
between A, p and h^^^ (and its derivatives). It can also be seen from (ll2l) that in 
the case = the Green's functions for A are those of Laplace's equation, rather than 
those of Helmholtz's equation. In this case A becomes long-ranged. In what follows we 
will therefore consider the ISS and kq = cases separately from the more general case, 
as they clearly have different behaviour. 

We can now integrate Eq. ( I39l) using the Green's function of Laplace's equation. 



The massless mode H^qq is then given by 



hS! = 2G{m + nXlfU 



47r(l + XI 



V 



M^Xl 
8vr(l + XI) 



V{h 



(m)2 
'■Oi 



i(m)Am) 



where the Newtonian potential U is 



U 



p{x') 



\x — X' 



■d^x' 



(44) 



(45) 
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and we have defined the additional functional 



V(0) = / , ' \, dV. (46) 
J \x — x'\ 

If M is large, we are then left with the long-ranged part of this metric component being 
given by 

h^Qfj=2G{m + nX^fU. (47) 

Here, and throughout, the symbol = will be taken to mean 'equal up to exponentially 
suppressed terms'. Such suppression is guaranteed as long as is much smaller than 
the length scale over which observations are being made, and can be arranged by choice 
of the parameters kq, ki, K2, / and Xq. 

Let us now consider the massive modes. These are different, depending on the 
theory in question. 

ISS theory 

For the ISS theory we can write the solution to ( HOil and ( H2|l as 

(^) _ 8G{m + n(r){m~n) (11 - 9a) / 

hoo ^ ^^^(^) + 247r(l + a)^"H^' 

(1-30-) /,M,{m)\ , (l-4a) 



247r(l + (T) V ""'^ ^^'V 6nl 
where we have defined the new functional 



H-M (l^h^) + ^^Wm (/4"'=) , (48) 



Wc(^) = / j^^^e-'"-' 'd'x', (49) 

J \x — X \ 

and suppressed the exponentially increasing mode. For sufficiently large M, the long- 
ranged part of this massive mode is zero: 

/^io^=0. (50) 

Kq = theory 

When kq = the solution to equations PUI) and fj42l) is 

h^m^ = ^(^ + nX',)im - n)WM (p) - j 



2\ 



^ (1-^0 



47r(l + Xi) 



(4(l-Xg)/., + (l-3Xg)/.2) / M2\ .... 

+ y^M[K, ), (51) 
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_ 8Gim + nX',){m-n) ^^ _ ((3 - Xg)Ki + (1 - Xg)/.^) ,, / ^2 



37r(l+X2)(2/€i + /€2) 



Oi,j Oj,i 



4(8., + (3-Xo^)._2)^. (.)2X ^52) 



37r/2X2 

and where we have used 3M^>Vm(^) = 12tcA — 32TcG{m + nXl){m — n)WM{p)- From 
the above it can be seen that even for large M the massive modes are not short-ranged: 

h^^'^=-—{m + nXl){m-n)U. (53) 
3 

Other theories 

For all other theories the solution to equations (HOl) and (H2l) is 
= + nXo^)(- - -)Wm (p) - ^^^^1^^ 



+ 



4vr(l + X2) 
(Xo^«:o - 8(1 - Xo^)/€i - 2(1 - 3X2)«:2) 



47r/2X 



Wm (/^r^) , (54) 



where 



_ 2GXgM^/^(m + nXg)(m-n) 
(X2ft:o + 12«:i + 6/^2) (1 + X^) ''^^ 

(X2(l - X2)/€o + 8(3 - X2)«:i + 8(1 - X^)/^^) 



47r(X2/s:o + 12fi:i + 6^2) (1 + X^) 
(Xo2(l - Xl)K, + 8(1 - Xl)K, + 8^:2),^, M, M 



47r(l + X2)(X2«:o + 12«:i + Qk,) '^^^ I'^o^.^' ''o,^ 
M^(8/.i + (3-Xg)K2) / M2A 

7r(l + X2)(Xo% + 12«:i + 6«:2)^V^°^ ^ ^ ^ 

and where we have defined 

2 ^ 2/to(l + Xg)(Xg/to + 8/€i+4/t2) 

/2(X2«o + 12«:i + 6/€2) ' ^ ' 

and used 

?Mi±^»l^ = 4G(m+„A-)(m-„)W„(rt + 5^^^m^i^W„(^).(57) 

Now, if both M and X are large enough then there are no long-ranged terms in the 
massive mode: 

C^=0. (58) 



The Parameterised Post- Newtonian Limit of Bimetric Theories of Gravity 



13 



5. Post-Newtonian Perturbations to 0(2) and 0(3) 
5.1. The gij and qij terms, to 0(2) 

Let us now consider the Qij and qij modes to 0(2), which we call 'post-Newtonian tensor 
modes', as they are tensors under transformations of the spatial coordinate system. 
From now on we will not consider further the contribution of the /iq™'^ ~ 0(1) terms to 
the metric perturbations. If they do not have an effect at the Newtonian level, then it 
seems unlikely they will be required at higher orders. 
Now, if we apply the gauge condition 

h^'^^,, - = (59) 

then the massless i-j equation, to 0(2), becomes 

V'/i^f = -87rG(m + nX^fdi.p, (60) 
which has the solution 

h[f = 2G{m + nXlf5ijU. (61) 

The massive modes, however, are once again dependent on the theory that is being 
considered. Now, the i-j equations, to 0(2), are 

V'/il^^ - M^/iJ^) + 87rG(m + nXl){m - n)5ijp 

2(2/^1 + K2) 2(2/^1 + /.2)(1 + Xg) ^ ^ 

" iXiKo + 8K^ + AK,)''' PXi ^^^^ 

where we have used the Bianchi identities to 0(2), Equation f H5]) . and where A is again 
given by the solution of Equation P2|) . 

ISS theory 



In the ISS theory the solution to f l42|) and fl62|) is given by 

(m) 40 4G 
% = -i7Tj:^{'m+na){m-n)WM{p,ij)+—{m+na){m-n)6ijWM{p)-{Qi) 



The term involving WA/(p,jj) is not in keeping with the usual PPN philosophy of 
avoiding terms with derivatives of p involved. Normally, one could perform a gauge 
transformation to remove such terms. Here, the massive modes are gauge invariant, 
but in subsequent sections we will find it possible to remove these terms from the 
combination of metrics that couples to the matter fields. For large M we have 

h^^=0. (64) 
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Kq = theory 



When Ko = the solution to and ([62]) is 



2G AG 
^iT^ = — ("^ + - n)^ijU + — (m + nX^){m - n)6ijmi{p) 

G 



m + nX^){m-n)WM{UAj). (65) 



371 

There is again an unsightly term containing derivatives of its integrand, which we will 
gauge transform away later on when considering the combined metric Qij. The long- 
range component of h[^^ is now given by 

h^^ = '^(m + nX^)im - n)5ijU. (66) 

Other theories 

For all other theories the solutions to ( H2l) and ( l62l) are given by 

AG 2G 

ht^^ = — (m + nXl){m - n)6ijWM{p) + -^{m + nXl){m - n)6ijWN{p) 



where A is given by ( l55l) . Once again it can be seen that in the general case there are 
two different massive modes, with masses M and N . Once again, there is also a term 
containing derivatives. If both M and are large enough then all of the terms in this 
massive mode are suppressed: 

/i!7^=0. (68) 
5.2. The goi and qoi terms, to 0(3) 

Now consider the g^i and go« modes to 0(3), which we call 'post-Newtonian vector 
modes', due to their properties under spatial coordinate transformations. If we again 
apply the gauge condition fl59l) . together with the new condition 

h^''\, - \h^'^\,o = -^C!o, (69) 
then we can write the massless combination of the field equations ([5]) and ([6]) as 

V^^o) = j{m + nX^fV^V, - W,) + 167rG(m + nX^^pv,, (70) 

where we have used U^i = |V^(H^i — Vi), as well as the previously found solution for 
h^Q . It can then be seen that equation (170|) has the solution 

= -^(m + nX'.rV, - f + nX'.rW,, (71) 
where Vi and Wi are the usual vector post-Newtonian potentials, as defined in 
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If we now use the 0(2) Bianchi identities, (143 p . as well as the Bianchi identities to 
0(3), 

2(X2fi:o+8fi:i+4/s:2)4S = (^o'^o+4/€i+2/€2)/iS + (X2/to+12/s:i+6K2)/iK, (72) 
then we can write the relevant combination of field equations from ([5]) and ([6]) as 

(^AqKo + OKi + 4K2j 

As above, the form of the massive modes in g^i depend on the theory being considered. 
ISS theory 

The solution to Equation (!73|) in the ISS theory is given by 

hr = , \ '-WMip,oi)-4:G{m+na){m-n)WM{pVi).{74:) 

3(i + crj 

As with the 0(2) components, there is still only one mass, M. If M is large enough 
then 

h^-^=0. (75) 

Ko = theory 

When Ko = the solution to Equation (!73|) is given by 

6(1 + Xq)[2Ki + K2) 

_ GfXl{m + nXl){m - n) 

6(l + X2)(2/€i + /€2) 

Possible long-ranged forces are now given by 

^(^)^ 0/^Xg(m + nXl){m-n) ^^ 

- Q{l + Xl){2^, + n,) ^^^^ 

Other theories 

For all other theories we have that the solution to (|73ll is given by 

^(^)_ 2^o(l + Xg) ^ 40(m + nXg)(m-n) ^^, 

% - 3/2Ar2M2 3M^ '^M(p,o.j 

-4G(m + nX2)(m-n)>VA/(pi;,). (78) 
To 0(3) we have no long-ranged terms in h^f^\ if M and are both sufficiently large: 
h^^^=0. (79) 



>Vm(p,o.). (76) 
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6.1. The gff and q^^' terms, to 0(4) 

Using known solutions and identities we can write the massless combination of the field 
equations ([5]) and (E]) as 



(0) 



2^(0) 



00 



-8TTG{m + uX^Yp 



l + U + 2v^ + 3 



MGXlim + nXl){m — n)ph, 



h. 



(0)2 
00,i 



-v-2/'/,(»")2 
^0 V''00,i 



Am) Am)-. 
"'00,ij"'ij ) 



+ 



2{2Ki + K2) 



00 



M2/2 



'^00,i^>« 



(1 + X^o) 

This equation can be integrated to give 



(m)2 



(80) 



(0) 
00 



2{m + nXlfG 



U 



(m + nX2)2 im + nXlf 



-GU^ + 2$i + 2- 



G$2 + *3 + 3$4 



2GX2(m + nX2)(m 



XoMV(/i^")' 



(1 + xi 

(m) 
00,i 



n],., . (m) ^ 
'00 / 



-V(p/i, 



v(/.K,/.!f; 



4vr(l + X2) 
(X^Ko + 12^1 + 6K2) 



+ 



XqKo 

(2/ti + K2 
27rM2/2 



2 V(4o^') 



■v(4S^.) 



167r/2 



V(A2 



167r(l + X2) 



^1) 



where the h^^ and A to 0(2) are taken to be given by the expressions found in 



''00 ' "'ij 



previous sections. The long-ranged modes in h^J are then given, for our various different 
theories by the following: 



ISS theory 



hS!=2im + n(rfG 



U - 



[m + na) 

a + a) 



AixaG'^l'^irn + naYim 
3(l + a)2 



n] 



1 + cr 



-G^2 + $3 



3$4 

(82) 
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K,o — theory 



9(1 + X2)M2 ' ' 3M2(1 + X2) 

+ 3M2(1 + X2) nP,^M{p),). (83) 



Other theories 



(84) 



In deriving these expressions we have used |Vf/p = ^V^C/^ — V^$2, as well as 

A^iB^i ^-V^{AB) + 2'KaA + 2'KpB--{M'^ + N^)AB (85) 

AiJB,^J = ^V^A^iB^i) + 2TraiA^i + 27rp,i5, - ^(M^ + N^)V\AB) 

+ ^{M^ + N^AB - tt{M^ + N^){Aa + Bp) (86) 

where A — Wm{p) and B — >VAr((T). Further manipulation has also been done by noting 
that, for example, V(p,j>VM(/o),i)=47rV(p^). These manipulations are performed so that 
the potentials are written without derivatives acting on the matter fields, as is usual in 
the PPN formalism. 

6.2. The c/^^^ and gj^^ terms, to 0(4) 

The equations involved in finding g^^^ and q^^^ to 0(4) are very lengthy, and we therefore 
choose to present them in Appendix B. Here we only state the form of the equations, 
and give their long-ranged components when M and N are large. We find the equation 



(m) 

2/,M _ A/r2^M 



for Hqq looks like 



V'h'^' - M%'^> = .... (87) 

The Green's function for this equation is therefore that of Helmholtz's equation, with 
the same mass term, M, as at 0(2). The right-hand side of this equation has terms 
containing A to 0(4). This quantity is given in Appendix B by the solution to an 
equation of the form 

{XIkq + 12ki + 6k2)VM - M'^XIkqA = .... (88) 

The Green's function for this equation can also be seen to be the same as the 
corresponding 0(2) case. The long-ranged components of the massive modes of Hqq 
to 0(4) are given below. 
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ISS theory 



/,M = ^L!i!^(^ + na^im - nf . (89) 



/to = theory 



In this case the massive modes /lo™"* and /i^-J^^ to 0(2) both have long-ranged 
components, and so the full expressions do not simplify greatly when the exponentially 
suppressed modes are neglected. The reader is therefore referred to equations fll35p and 
(IT39I) in Appendix B, where to 0(4) is given explicitly. 

Other theories 

4o^=0. (90) 

To find these expressions we have made use of fISS]) and fl5^ . as well as the following 
relations: 

Wm(>V^(X)) = (W^(X) - Wm{X)) (91) 

^m{V{X))~V{X). (92) 

7. Perturbations Coupled to Matter Fields 

We have so far calculated the form of the massless and massive combinations of the 
perturbations to g^i, and q^y. We now want to know what these tell us about the 
perturbations to the metric g^u. Using previous relations we can write g^y as 

gfiu = Vfiu + (93) 



where 



^ _ (m + nXg) Xg(m-n) 
^'^ (1 + X2) (1 + X2) ^ ' 



and fj^y = rj^y + rj^^. It is now straightforward to write down the components of h^j^y. 
To proceed in doing this it is useful to define two new constants: 

_ Gjm + nXlf 
= (1 + XI) 

Qn - ^^^^^ [3(m + nXlf - Xl{m - r^f] . (96) 

The reason for doing this is that although the constant G appears in the gravitational 
action ([2]) in the same way as it does in general relativity, it is not necessarily the value 
of Newton's constant that is determined by Cavendish type experiments. For the = 
theories this value of Newton's constant is given by (when M is large), while for the 
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ISS theory and all other theories it is given hj (when M is large for the ISS theory, 
and when both M and are large for all other theories). 

It is also useful at this point to start using infinitesimal gauge transformations 
to put our results in the 'standard post-Newtonian gauge'. In the preceding section 
we used the gauge specified by conditions (|59l) and ( |69ll . This has been a convenient 
choice, and has allowed integration of the field equations to post-Newtonian accuracy. 
It has, however, resulted in the metrics in question having off-diagonal components in 
the spatial part of the metric, as well as non-zero 0(1) terms in the 0-i components of 
the metric. These can be removed by making an infinitesimal coordinate transformation 
of the form x'^ ^ x'^ + The metric that couples to matter, g^^iy, is then transformed 
in such a way thatl^ 



where covariant derivatives are with respect to gf^i^, and where indices have been lowered 
with g^i,. By making coordinate transformations such that ~ 0(1) or 0(3) and 
C,i ~ 0(2) the metric perturbations then transform as 



Transformations of this kind will be used in what follows to remove all contributions 
from the 0(1) massive modes in hoi, as well as to diagonalise hij, and to remove the 
potentials found above that contain derivatives on the functions in their integrands. 

As in the previous section, many of the expressions involved here are quite lengthy. 
We will therefore once again present the full equation in Appendix C, and quote here only 
the terms that are long-ranged when M and are large (after the gauge transformations 
discussed above have been performecill). To 0(4) in hoo, 0(3) in hoi and 0(2) in hij we 
find the results below. 



(97) 



hij }■ hij 2^(^jj^ 
hoi — > hoi — ^Q^i — ^ifi. 



(98) 
(99) 



ISS theory 



hoo—2GN U 




+ 2$i + 




$2 + $3 + 3$4 



3 V(P^) + 



TrH^Glajm - n)^ 
18(1 + af{m + na) 



P 



2 



(100) 



■Oi — 



1 2 



(101) 



hij—2Gi\f6ijU. 



(102) 



% This is true as long the matter does not couple to the massive combination of metrics only, i.e. as 
long as m ^ —uXq. 

+ These transformations are given explicitly in Appendix C. 
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Kq = theory 



hoo=2gNU + 0(4) 



3(m + nX2)2 



^ 3(m + nXg)^ 

2 (3(m + nX2)2 -Xg^fm 



■ n 



'(3(m + nX2)2-X2| 



(103) 

(104) 
(105) 



Other theories 



^00 = 

hoi= 
hij= 



-2Gn 
7Gn 



U - 



G 



N 



[m + uXq 
Gn 



-f/2 + 2$i 



2G 



N 



[m + nXn 



3$4 



2 

-2G i\f6ijU. 



(106) 

(107) 
(108) 



In the Kq = theory we have not written the 0(4) part of /iqo explicitly, because, as 
mentioned before, it contains a large number of long-ranged potentials. Appendices B 
and C contain enough information for the reader to calculate these terms explicitly, if 
they are required. 



8. Post-Newtonian Parameters, and Observational Constraints 

It is convenient at this point to choose units so that m + nX^ = 1. The unperturbed 
metric fj^u then takes its usual form. Using the results of the previous section we can 
now obtain the PPN parameters of the theories we are considering, and determine the 
extent to which they can be observationally constrained. 

ISS theory 

It can be seen by comparing ffTOOj) . (fTOTj) and (11021) with (EID, (122]) and ([23]) that, 
when M is large, all PPN parameters take the same value as in General Relativity. We 
do, however, have two new potentials in the /iqo term at 0(4). These are 

3(1 + <t) ^ 18(l + af ' 

No such terms appear in the usual PPN prescription. The second of these is particularly 
unfamiliar as it involves p directly, and not p integrated over any volume. Such a term, 
however, can be expected to be small in comparison to V(p^) so long as / ^ rf, where d is 
the length scale over which the observational phenomena are being measured. Outside 
of any matter distribution, of course, this term will vanish entirely. 
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In the absence of any observational constraints on these new potentials, the only 
constraint that can currently be placed on this particular theory is therefore Md ^ 1, 
or equivalently 

4(1 + a) »^. (110) 
No other constraints are available from standard PPN formalism. 

kq = theory 

In the Kg = theory it can be seen by comparing fll03p . fll04p and fllOSp with fl^ . 
(!22|) and (123|) that, among other deviations from general relativity, we have 

(3+xg(i-n(i-xg)r 

{3 - Xiil - n{l - Xi)y 

As in the ISS theory, the = theory must also obey Md ^ 1 in order for the massive 

modes to be exponentially suppressed. For this theory, this corresponds to 

8(1 + X2), , /2 

^ ^\2k, + K2):>-j.. (112) 



l-^^^Ol^-'H^-^o;; ) ^^^^ 



Unlike the ISS theory, however, we can also use the PPN constraint on 7 from 
observations of the Cassini spacecraft to constrain the theory. To be compatible with 
observations 7 must satisfy equation ( 12^ . We then have the constraint 

X^{1 - n{l - X^))^ = (3.1 ± 3.3) X 10"^ (113) 

While tight, it is worth noting that this constraint can always be evaded by an 
appropriate choice of Xq. Unless Xq is known from other considerations (cosmology, 
for example) it is not possible to constrain the theory with this information only. It is 
also worth noting that this is a constraint on n only, and not the parameters determining 
the interactions in the theory (i.e. ki and ^2)- Further constraints are also available on 
this theory from the ^00 metric component at 0(4). 

Other theories 

Comparing ffT06D . ffTOTj) and ffT08|) with ([2l]), ([22]) and ([23]) shows, as long as both 
Md 3> 1 and Nd ^ 1, that all of the PPN parameters in this general class of theories 
are the same as in general relativity. The only constraints that need to be satisfied by 
these theories are therefore 

^^^^f'\ x',Ko + 8k, + Ak,) » ^ (114) 

and 

2"°^' + ^°)(X^o + 12., + 6.2) ^'''^ 
If these conditions are met then all the post-Newtonian predictions of these theories are 
indistinguishable from those of general relativity. This shows the importance of the kq 
term in the interaction Lagrangian. 
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9. Discussion 

We have investigated liere tlie weak field limit of a general class of bimetric theories, 
generated from the gravitational action ([2]), and coupled to matter via a linear 
combination of the two metrics involved. We have calculated the post-Newtonian limit 
of these theories, up to 0(4) in qqq, 0(3) in goi and 0(2) in gij, for a perfect fluid matter 
content distributed without any symmetries. We find that, to all orders considered, 
there is a natural decomposition of the perturbations to the two metrics into modes we 
have called 'massive' and 'massless'. These modes are defined through the combinations 

/i!.'^^ = V-% and /iW = V + V, (116) 

^0 

where h^j^^ and h^i, denote the perturbations to the two metrics involved, and Xq is a 
constant associated with the unperturbed background space-times. 

We find that, in general, the 0(1) contributions to the massive modes /iq™'' are non- 
zero, and cannot be made to vanish by a coordinate transformation. These quantities 
satisfy a Helmholtz equation with mass M (as given by Eq. ( l35i) ). and while the 0-i 
component of the combined metric that couples to matter can always be set to zero 
by a gauge transformation (as long as the matter does not couple exclusively to the 
massive modes), these 0(1) terms can in principle make a non-negligible contribution 
to the other components of the metric. We therefore include the effects of these terms 
while calculating the 0-0 component of the combined metric to 0(2), but neglect them 
to higher orders, as we see no reason to expect them to contribute substantially to 
post-Newtonian accuracy if they are negligible at the Newtonian level. 

To Newtonian order, we find that in general the relevant metric perturbations 
have massless modes (that are always long-ranged), as well as massive modes with two 
different masses: M and (where is given by Eq. (I56l)). We find two special cases 
in which the second mass, A^, does not appear. These are the ISS theory, as specified 
in Section El and theories with = 0, in which the interaction term yj—q{q~^)^^g^iv is 
missing from the action. In the ISS theory this occurs because the V/ig^'* = ... equation 
can be decoupled from the h^^^ terms without having to solve any other differential 
equations. In all other cases this is not possible, and a second equation of the form 
V^A — N'^A = ... must be solved, before one can integrate the /iq™'* equation. In the 
second anomalous case, when — )■ it is found that the mass parameter N ^ 0. 
Theories of this kind therefore have extra massless modes, beyond the usual Newtonian 
term that occurs from integrating the Vh'^Q = ... equation. 

We then proceed to calculate the metric perturbations to post-Newtonian accuracy. 
These terms all involve the same two masses, M and A^, that are required to Newtonian 
order (except for the ISS and kq = cases, which only require M). By expressing these 
perturbations in the standard post-Newtonian gauge we find that the predictions of 
these theories are, in general, very similar to those of general relativity. In particular, if 
M and A^ are both large enough then they can evade all preferred frame experiments, 
and are indistinguishable from general relativity by any experiments that constrain the 
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PPN parameters /3 and 7 (including the stringent tests available from Shapiro time-delay 
observations [39], and lunar laser ranging measurements [IQ]). The ISS theory also does 
a remarkably good job of satisfying weak field tests of gravity when M is large, with 
PPN parameters that are again the same as in general relativity. In this case, however, 
we do find a couple of anomalous terms appearing at 0(4) in the 0-0 component of 
the metric that couples to matter. These are given in Eq. (11091) . and could potentially 
be used to observationally distinguish between general relativity and ISS theory. The 
Ko = theory does not fare quite so well, due to its extra long-ranged modes. In this 
case observations of 7 can be used to place constraints on a particular combination of 
the constant Xq, and the mixing angle between the two fundamental spin-2 fields that 
couples to matter, even when M is large. 

As well as the large mass limit, one can also consider the small mass limit. In this 
case a viable Newtonian term can again appear in the 0-0 component of the metric, 
and we must then consider the problem of the vDVZ discontinuity that was discussed 
in Section [H To do this we note that as M — )■ and — )■ we get Wm(p) — ^ U* and 
Wn{p) — ^ U\ respectively (we have used the superscripts * and f here, even though 
f/* = f/t = so that we can keep track of where the potential has come from). 
Neglecting contributions from h^^^ at 0(1), the full expressions in Appendix C then 
give, for the most general theories, 

t^^-nC,,U I + 20^.Xg(l-n(l + Xg))^[/t ^^^^^ 



h,, = 2Gr,U + ^GNXl{l-n{l+Xl)fU* ^ 20^Xg(l - n(l + Xg))^[/t _ ^^^^^ 

The vDVZ discontinuity can be seen to be manifest in the U* terms of these two 
equations. If these were the only terms in these expressions then we would indeed 
get the familiar 7 = 1/2 that is anticipated from massive gravity. Here, however, we 
have two additional terms in both Iiqq and hij - one from the terms that are long-ranged 
even when M and N are large, and another from the VVAr(p) terms. The cumulative 
effect of all three of these contributions added together shows that there is no vDVZ 
discontinuity in the 7 parameter of the most general class of these theories, as in this 
case we end up with h^ = Sijhoo, and hence 7 = 1. The same is true in the special 
case in which = 0, as should be expected, as these modes already correspond to the 
X — )■ limit of the more general theories. The ISS theory, however, does not have 7 = 1 
as M — )■ 0. Instead we find 

3 + 2(l-n(l+.y. (119) 
' 3 + 4(1 -n(l + a))V ^ ' 

This is essentially because the WnIp) terms that cancel the discontinuity from the 
Wm(p) terms are absent in the ISS theory. 

It should be noted, however, that although the general relativistic value of 7 is 
recovered in the most general theories as M and X — )• 0, this does not mean that all 
PPN parameters will be. For example, for the parameters ai and 0:2 that parameterise 
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deviations from general relativity in the 0-i components of the field equations, we find 
that while ai — )■ (as in general relativity), we also have 

(1 - n(l + Xi)m 

^ -{2-„{i + xs)y-xi 

In terms of this expression, equation ( 1251) then gives the constraint 

Xo2(l-n(l + Xo2))2< 1.2x10-^. (121) 

This is the same combination of quantities that was constrained in the large M limit 
of the kq = theory using the parameter 7, but the constraint here is from an entirely 
different gravitational effect, and is even tighter. On should also note that in the limit 
M — the 0(1) terms in /iq™'^ are no longer short-ranged. These terms are, however, 
set by initial conditions, rather than by the matter content of the space-time. One may 
then be able to limit their magnitude by applying suitable observational constraints, but 
this will place constraints on the initial conditions only, and not on the gravitational 
theory itself. 

Finally, we note that when the matter fields are coupled to the massless combination 
of metrics only (so that m = n) then many of the constraints found above are 
automatically satisfied. These observations can therefore be though of as constraining 
the amount of massive mode that is allowed to be mixed with the massless modes. 



Acknowledgements 

TC acknowledges the support of Jesus College, Oxford and the BIPAC. MB would like 
to thank the Oxford Astrophysics department for hospitality during a sabbatical year. 
MB was partially supported by the JS Guggenheim Memorial Foundation, Fondecyt 
(Chile) Grant # 1100282 and Alma-Conicyt Grant # 31080001. 



The Parameterised Post-Newtonian Limit of Bimetric Theories of Gravity 



25 



References 

[1] C. Brans, and R. H. Dicke. 'Mach's principle and a relativistic theory of gravitation'. Phys. Rev. 
124, 925 (1961). 

[2] C. M. Will, and K. Nordtvcdt. 'Conservation laws and preferred frames in relativistie gravity. I. 

Preferred-frame theories and an extended PPN formalism'. Astrophys. J. 177, 757 (1972). 
[3] C. J. Isham, A. Salam, and J. A. Strathdee. 'F-dominance of gravity'. Phys. Rev. D 3, 867 

(1971). 

[4] K. Nordtvedt. 'Equivalence principle for massive bodies. II. Theory'. Phys. Rev. 169, 1017 
(1968). 

[5] C. M. Will. 'Theoretical frameworks for testing relativistic gravity. II. Parameterized post- 
Newtonian hydrodynamics and the Nordtvedt effect'. Astrophys. J. 163, 611 (1971). 

[6] C. M. Will. 'Theory and experiment in gravitational physics, revised edition'. Cambridge 
University Press, Cambridge (1993). 

[7] N. Rosen. 'A bi-mctric theory of gravitation'. Gen. Rel. Grav. 4, 435 (1973). 

[8] D. L. Lee. C. M. Caves. W.-T. Ni, and C. M. Will. 'Theoretical frameworks for testing relativistic 
gravity. V. Post-Newtonian limit of Rosen's theory'. Astrophys. J. 206, 555 (1976). 

[9] P. Rastall. 'A theory of gravity'. Can. J. Phys. 54, 66 (1976). 
[10] P. Rastall. 'A note on a theory of gravity'. Can. J. Phys. 55, 38 (1977). 

[11] A. P. Lightman, and D. L. Lee. 'New two metric theory of gravity with prior geometry'. Phys. 
Rev. D 8, 3293 (1973). 

[12] C. M. Will, and D. M. Eardley. 'Dipole gravitational radiation in Rosen's theory of gravity: 
Observable effects in the binary system PSR 1913+16'. Astrophys. J. 212, L91 (1977). 

[13] T. Damour, I. I. Kogan, and A. Papazoglou. 'Non- linear bigravity and cosmic acceleration'. 
Phys. Rev. D 66, 104025 (2002). 

[14] M. Bahados. 'Eddington-Born-Infeld action for dark matter and dark energy'. Phys. Rev. D 77, 
123534 (2008). 

[15] M. Banados, P. G. Ferreira, and C. Skordis. 'Eddington-Born-Infeld gravity and the large scale 

structure of the Universe'. Phys. Rev. D 79, 063511 (2009). 
[16] M. Bafiados, A. Gomberoff, D. C. Rodrigues, and C. Skordis. 'A note on bigravity and dark 

matter'. Phys. Rev. D 79, 063515 (2009). 
[17] C. Skordis. 'Eddington-Born-Infeld theory and the dark sector'. Nucl. Phys. Proc. Suppl. 194, 

338 (2009). 

[18] M. Milgrom. 'Bimetric MONO gravity'. Phys. Rev. D 80, 123536 (2009). 

[19] T. Clifton, and T. Zlosnik. 'FRW cosmology in Milgrom's bimetric theory of gravity'. Phys. 

Rev. D 81 103525 (2010). 
[20] N. Boulanger, T. Damour, L. Gualtieri, and M. Henneaux. 'No consistent cross-interactions for 

a collection of massless spin-2 fields'. Physics AUG 10, 94 (2000). 
[21] T. Damour and 1. 1. Kogan. 'Effective Lagrangians and universality classes of nonlinear bigravity'. 

Phys. Rev. D 66, 104024 (2002). 
[22] D. Bias, C. Deffayet, and J. Garriga. 'Global structure of bigravity solutions'. Glass. Quant. 

Grav. 23, 1697 (2006). 

[23] D. Bias, C. Deffayet, and J. Garriga. 'Bigravity and Lorentz-violating Massive Gravity'. Phys. 

Rev. D 76, 104036 (2007). 
[24] D. Bias. 'Aspects of Infrared Modifications of Gravity'. PhD thesis (2008). 
[25] Z. Berezhiani, D. Comelli, F. Nesti, and L. Pilo. 'Spontaneous Lorentz breaking and massive 

gravity'. Phys. Rev. Lett. 99 131101 (2007). 
[26] Z. Berezhiani, D. Comelli, F. Nesti, and L. Pilo. 'Exact Spherically Symmetric Solutions in 

Massive Gravity'. JHEP 07, 130 (2008). 
[27] Z. Berezhiani, F. Nesti, L. Pilo, and N. Rossi. 'Gravity Modification with Yukawa-type Potential: 

Dark Matter and Mirror Gravity'. JHEP 07 083 (2009). 



The Parameterised Post- Newtonian Limit of Bimetric Theories of Gravity 26 

[28] D. Bias, D. Comelli, F. Nesti, and L. Pilo. 'Lorentz Breaking Massive Gravity in Curved Space'. 

Phys. Rev. D 80, 044025 (2009). 
[29] H. van Dam and M. Veltman. 'Massive and massless Yang-Mills and gravitational fields'. Nud. 

Phys. B22, 397 (1970). 

[30] V. I. Zakharov. 'Linearized gravitation theory and the graviton mass'. JETP Lett. 12, 312 
(1970). 

[31] D. G. Boulware, and S. Deser. 'Can gravitation have a finite range?'. Phys. Rev. D 6, 3368 
(1972). 

[32] W. Pauli, and M. Fierz. 'On relativistic wave equations for particles of arbitrary spin in an 

electromagnetic field'. Proc. Roy. Soc. Land. A73, 211 (1939). 
[33] S. L. Dubovsky. 'Phases of massive gravity'. (arXiv: hep-th /0409 1 24| 

[34] V. A. Rubakov, and P. G. Tinyakov. 'Infrared-modified gravities and massive gravitons'. 
arXiv:0802.4379 [hep-th]. 

[35] A. H. Chamseddine and V. Mukhanov. 'Higgs for Graviton: Simple and Elegant Solution'. 

larXiv:1002.387 7 [hep-th]. 
[36] A. I. Vainshtein. 'To the problem of non-vanishing gravitation mass'. Phys Lett. B 39, 393 

(1972). 

[37] N. Arkani-Hamed and H. Georgi and M. D. Schwartz. 'Effective field theory for massive gravitons 

and gravity in theory space'. Ann. Phys. 305, 96 (2003). 
[38] E. Babichev, C. Deffayet, and R. Ziour. 'Recovering General Relativity from massive gravity'. 

Phys. Rev. Lett. 103, 201102 (2009). 
[39] B. Bertotti, L. less and P. Tortora. 'A test of general relativity using radio links with the Cassini 

spacecraft'. Nature 425, 374 (2003). 
[40] J. G. Williams, S. G. Turyshev and D. H. Boggs. 'Progress in lunar laser ranging tests of 

relativistic gravity'. Phys. Rev. Lett. 93, 261101 (2004). 
[41] K. Nordtvedt. 'Probing gravity to the second post-Newtonian order and to one part in 10^ using 

the spin axis of the Sun'. Astrophys. J. 320, 871 (1987). 



The Parameterised Post- Newtonian Limit of Bimetric Theories of Gravity 



27 



Appendix A: Perturbed Metric Quantities 

The perturbed form of the metric that we will be using is specified in (|TOl) and (iTTl) . In 
order to solve ([5]) and IQ, however, we also need to know the perturbed Ricci tensor 
to 0(3) in its off-diagonal components, and 0(4) in all of its diagonal component^], 
as well as the relevant expressions for the perturbed metric determinants, and inverse 
metrics. The perturbed Ricci tensor components take on their usual lengthy form, and 
we will not reproduce them explicitly here. 

The perturbed metric determinants can be calculated using the equation det{g) = 
Qtr{in{g))^ To 0(2) in perturbations this gives 

det{g) = - (1 - /loo + ha + hl^) + 0(4) (122) 



c/e%)= -X«(l-^ + || + ||)+0(4) (123) 

-^0 -^0 ^0. 



det{g) 



(m + uXq^ + {ha - hoo + hli)m'^ 
+ {Xlhii + ha - XqHqo - hoo + 2hoihoi)mn 



+ {xlh,, - xlhoo + nl,)n^ 



m + nXo')' + 0(4), (124) 



and, if we ignoring 0(1) contributions to hoi and hoi-, then to higher order we have 

2 2 

hoo ha hoohii hf; h'fn 



det{g) = -(1-hoo + hu- hooh^ + -hi - ) + 0(6) (125) 



J +1 \ v^f ^ I hoohii h% h^j \ noR\ 

det[q) = -Xq I-— 2+— 2 - -^rr- + ^ - ^ + 0(6). (126) 

In order to compute the interaction terms in the field equations ([5]) and ([6]) we also 
need to know the inverse metric fluctuations {q~^)'^'' = {fj~^)^'^ + {6q~^)^'^ . These can 
be found from the definition {q~^)'^°'qau = ^u, which gives 

(5q'T = -irTiv-y'^aP + {v-T{v-'r\v-y^'^aP^^5 + 0{h^). (127) 
The components of {dq^^^ are then given by 



-uoo "oo , 



hoQ_ h^ 



+ -^ + 0(4) (128) 

^0 ^0 

-^^°'-^ + 0(3) (129) 



'7^ = -|i-^ + 0(4). (130) 
-^0 -^0 

Ignoring 0(1) contributions to hoi and hoi-, the higher order components of {Sq'^'^ are 

-l\00 ^00 ^00 I ri(a\ 

) = - w ~ V6 +0(6) (131) 

^0 ^0 

^'f = -|i + %|^ + 0(6). (132) 
-^0 -^0 

* The i-j components are usually only required to 0(2), but here we also need higher order terms in 
order to calculate the 0(4) part of the 0-0 components. 
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Using this expression, together with other known identities and lower order solutions, 
the field equations ([5]) and give 

V ILqq IVl ILqq 

2{1 + XI){2ki + K2) 



-A 



A{1 + XI){2k^ + K2) 



MH^Xl 



Am + h 



(m) 



-%-KG{m + nXQ){m — n)p 
8TcG{m + nX^){mX^ + n 



l + U + 2v^ + 3- 
P 



00,00 

2(2/ti + K2) JO) 



(1+^0^^ 



■ph^^'^ + ATiGmn{l + Xl)p (h^^'^ + h. 



(m) 

a 



(1 + ^0') 
(1 - ^0') 



''■00,ij"'ij 



"'00,1 



2(1 + X2 

[X^Ko + 8ki + 6K2] 

4^2 



2k 



'^'^00,j'''00,i 

(m) , (m) 
00,ii"'ii 



+ 



2{X^Ko + 4fi:i + 2^2) , (0), (m) 



''00 "-00 



r)l,{m)Am) i{m)i(m) 



3/l, 



(m)2 
00 



+ 2/i(™)/i!r) + /ir^-2/i. 



{m)2 



(m)2 



'■DO 



""■GO + 72 \^2 00 % + 



(4/ti + 3^:2) (m)2 



The solutions to this equation are given below. 



(133) 



ISS theory 



h 



00 



2G{m + ncr)(mcr + n) 



(1 + ^) 



Wm(p4?) 



-Gmn{l + a) Wm(p4?) + ^^Miph^T^) ^ ,^ ^ 
L J o7r(_i + cr 

+2G(m + na){m - n) \Wm{p) + ^/(pn) + 2WM(pt^') + 3Wa/ 
(WM(4?..,/.!r^) - 2Wm(/.S,<])) >VM(4nS 



1,00; 



47r/2 



a 



87r(l + a) 



>VM(/.sr)-2>VM(4;:!,/.sr^^ 



-2>vm(/.K/.!;5) + >v^^(^E^5S) 



(134) 
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"■00 



(1 + X2)(2«i + K2) ... , Wm(^,oo) ^ Wm(/i£AO 
— KVMi^j H H 



^-2G{m + nX2)(m - n) \Wm{p) + WmCpH) + 2>Vm(p^^') + 3>Vm(p)] 

2{m + nXl){mXl 



(0) 7 (m)^ 



47r(l + XI) 



ttPXI 



{iKi + 3^2) 



3(/ti + K2) 



277^2^2 



47r/2X2 



(my 
ii 



{4:Kl + 3^2) 



inPX^ 



(1 - XI) 

87r(l + X2) 



2WM(4S,/^!r^) 



2>VM(/.g/.Si;^) + >VM(/.SCi) 



(m) , (m)N 



Other theories 

(l+X2)(2/€i + /€2) 



(m) 
00 



27r/2Xo2 



. . (l + Xo2)(2/.i + /.2) 
Wm(^) + _,^o,o ^2 I^mIAoo) 



(135) 



ttMH^X^ 



WmTOo) , (2«i + /^2) 



47r 



+ 



27rM2/2x2 



+2^(771 + nX^){m - n) [Wm{p) + VVnip^) + 2Wm{pv^) + 3>Vm(p)] 



2^(771 + 7zX2)(7nX2 + 7l) (^), 



{1+Xi) "-VA^-OO ; 

-^77771(1 + X^) [>VM(p/i&^) + WM(p/ilr^) 

(>VM(Co!./y ) - 2>Vm(C<!)) + K,) 

47r(l + X2) 

{X^Ko + Aki + 2k2 



47r/2X2 







47r/2X2 



27r/2X2 



277^2^2 



477/2^2 



(1 - ^0') 



87r(l + X2) 

[X^Ko + 8Ki + 6^2] 

16^ 



>VM(<?)-2WM(<l,/.lf) 



3>VM(4?') + 2>VM(/^ff/^r) 



(136) 
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The expressions above should be understood to have H^^q , h^^\ /i-™'' and A given to 
0(2) by the lower order solutions found in previous sections. It now remains to find A 
to 0(4). To do this we use the field equations ([5]) and ([6]), as well as known identities 
and lower order solutions, to find 



(X^Ko + 12ki + 6k2)VM - M'^XIkqA 



"-00,00 



2(1 + X, 



0. 



(1 + ^0^; 



SnGMH^X^im + nX^){m - n) 



3mn(l - X^) + 
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(1 + ^0^: 



nfXl 
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00 
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(1 + ^0^) 
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(m) 



A\ + AV^A 



2/i, 



(m)2 



00,i 



00 ^ "-oo 
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Am)2 
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(137) 



The solutions to this equation are given below. 
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(m) 
''00,00 



-5/.gV2/.rM,,, 
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= theory 
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(139) 
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In the above we have made use of the Bianchi identities to 0(4), which are given by: 
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As long as m 7^ —nX^ we can always use a coordinate transformation to remove all 
contributions from the 0(1) massive modes to /loi, so that hoi — > 0. We will therefore 
consider here only the higher order contributions to h/^^- 



Newtonian perturbations to 0{2) 

To 0{2) in hoo we find the results below 
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^ (4(1 - Ag)^i + (1 - 3Ag)^2)(m - n) ^^,^^ / (^)2 



Other theories 



27ri2(l + A2) 
A2(m + nA2) 



(143) 



47r(l + A2)2 

M2A2(/// + //A2) 



87r(l + A2)^ 



-V 



V (430 - ^ (^S^S) 

(m-70(8/vi + (3- A2)k2) 



(430 



+ 



37ri2(l + A2) 



(430 



+ 



8G7vA2(m 



n] 



-y\^M (p) 



2GjvA2(m 



n] 



3(m + nAo2)2 
(l_X2)(m-n)Ao2 ^ 
47r(l + A2)2 



3(m + nAo2)2 



-W^ (p) 



Wm ( 43' 
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(m - n){X^Ko - 8(1 - X',)k, - 2(1-3X|)^ / 



47r/2(l + X2) 



(m - n)(Xg(l - X^)no + 8(3 - + 8(1 - Xg)/^2) 

127rM2/2(l + x2) ''''^ 



■ (144) 



(m - n)(X2(l - X2)ft:o + 8(1 - Xi)Ki + 8^2) 
^ 12ttMH^{1 + X^) 

Post-Newtonian perturbations to 0{2) and 0(3) 

To 0(2) in hij and 0(3) in h^i we find the results below. 

ISS theory 

3(1 + (T)(m + ncr)2 ' ' 3(m + n(7)2 

3{m-\-na)^ 



and 



hoi = ^V^i - -T^Wi t^Wm 

2 2 (m + nc)^ 

^ ' r>VM(p,oi) (147) 



3(1 + (7) (m + n(T)2 



The arrows in the expressions indicate the infinitesimal coordinate transformations 

. GNl'^crim — Tif , , ,^ 

ii = - ^.-.^ w ^ \2 >VM P ,i 149 
6(1 + (7)(m + ncr)^ 

GNl'^crim - nY , , 

6(1 + (j)(m + ncrj^ 



= theory 



? _ _ (3(m + nX2)2 + x2(m-n)2) ^ 
Aii,- - ^^N ^3^^ ^ ^^2)2 _ xl{m - n)2) 



"(3(m + nX2)2 - X2(i - n)^)'^^^'^^^^^ 

G.Xiim-r^r ^^^^^^^^ ^^3^^ 



TT (3(m + nX^y - X^(m - n)2) 
(3(m + nX2)2 + Xg(m-n)^) 
(3(m + nX^y - X^{m - n)^) 
AgNX^{m-nfSij 



' (3(m + nX2)2 - X2(m - n)2) 



>Vm(p) (152) 
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and 



^ AGXi{m + nXi){m-nr 7G(m + nXg)^ 



GPX^{m + nX^){m - nf _ G{m + nX^f 

6(1 + X2)2(2ki + K2) 2(1 + 

G'/2x4(m + nX2)(m - n)^ 



6(1 + X^y{2n, + K2) 



>Vm(p,oO (153) 



7G(m + nXg)3 _ G(m + nXg)^ 
2(1 + Xo^) 2{1 + Xi) 



^^^"("' + "^")("'^W.(p..), (154) 



(1 + ^0^) 

where the arrows indicate the coordinate transformations 

, ^ GNl^X^im-nf 

4(1 + X2)2(2/^i + /s:2)(3(m + nX^)^ - X^im - n)^) ^ ^ 



GNl^X^im 



n 



i2 



- 4(1 + XIY{2k, + K,){Z{m + nXlf - Xl{m - nf) " 3>Vm(p),o). (156) 
Other theories 



n 



i2 



n 



i2 



and 



f. 7Gn., Gn,,, 4GNX^{m-nf 



_ TGat,, Gn 4G^vX|(m-n)2 
where the arrows indicate the coordinate transformations 



We will not show the ftill expressions for hoo to 0(4) here as they are quite lengthy, and 
we do not feel that writing them out explicitly will add sufficient extra insight to justify 
their inclusion. The interested reader can calculate these quantities straightforwardly 
using the relevant expressions in Appendix B, and the coordinate transformations given 
above. 



